Typeset with jpsj2.cls <ver.l.l.l> 



Possibility of Horizontal Line Node in KFe 2 As 2 
Probed by Muon Spin Rotation 
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The anisotropy and temperature dependence of the magnetic field penetration depth in the 
iron-based superconductor KFe2As2 have been measured down to 20 mK by means of muon 
spin rotation. We have observed that the flux-line lattice forms triangular symmetry when 
we applied the field parallel to c-axis, suggesting the interaction between vortices should be 
isotropic. The temperature dependence of the penetration depth observd with an applied field 
parallel and perpendicular to c-axis are different. These results can be accounted for by a 
superconducting gap function with a horizontal line node in the basal plane. 
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The discovery of superconductivity at T c = 26 K in 
LaFeAsOi-^Fj, 1 and the improvement of T c above 50 K 
in other iron pniceides 2 ' 3 have stimulated a tremendous 
research effort to elucidate the superconducting pairing 
mechanism and symmetry in iron pnictide superconduc- 
tors. In oxygen-free system of AFe2As2 (A — Ca, Sr, Ba, 
and Eu), 4 ' 5 the emergence of superconductivity seems 
to be associated with the suppression of magnetic or- 
dering, which is accomplished either via application of 
pressure 5-7 or by chemical substitutions. 4, 8-11 In the op- 
timally doped Bai_ :E K x Fe2As2 (x ~ 0.4) with T c = 38 K, 
the Fermi surface consists of well-separated hole and 
electron sheets, which seem to exhibit good interband 
nesting, as found in many iron pnictide superconduc- 
tors. For example, the angle resolved photoemission 
spectroscopy (ARPES), 12 muon spin rotation (^SR), 13 
and NMR 14 have suggested the multiple full gap with 
sign change between these two different bands, namely 
s±-wave. Meanwhile in KFe2As2, which shows super- 
conductivity at T c — 3.5K, de Haas- van Alphen experi- 
ment and the accompanying band structure calculations 
carried out by Terashima et al. 15 have shown that the 
electron sheets centered at the X point of a Brillouin 
zone are replaced with small hole tubes, therefore, the 
interband nesting observed in optimally doped sample 
is not there anymore. This unusual electronic proper- 
ties was also confirmed by ARPES 16 and recent neutron 
measurement. 17 Other studies, i.e., NQR and specific 
heat, 18 thermal conductivity, 19 NMR 20 and penetration 
depth, 21 have suggested that the superconducting gap in 
KFe2As2 has line node in contrast to the nodeless gap 
suggested in Ban.6Kn.4Fe2As2. On the other hand, an 
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isotropic triangular flux-line lattice (FLL) was observed 
by the small angle neutron scattering (SANS) without 
FLL transitions up to H C 2 with the condition of H par- 
allel to c-axis. 22 This result suggests that an interac- 
tion between vortices is isotropic. According to the the- 
oretical approach, it is suggested that horizontal node 
(node in afe-plane) appears in the hole band at T point 
of a Brillouin zone in KFe2As2, 23 while some closed the 
nodal loops at the electron Fermi surface was suggested 
in relative superconductor of BaFe2(Asi_ l; P 2 ;)2. 24 Al- 
though various possiblities of line node in KFe2As2 are 
suggested, the reported experimental data are observed 
in polycyristalline samples or in single crystals only with 
the condition of an applied field H parallel to c-axis. In 
order to elucidate the superconducting gap structure in 
KFe2As2, we have observed the temperature dependence 
of penetration depth for both H || c and H _L c-axis 
down to 20 mK. 

In this Letter, we report the demonstration of the hor- 
izontal nodal superconductivity in KFe2As2 by means of 
/iSR measurements. Our findings reveal that the tem- 
perature dependence of penetration depth measured only 
with one field direction, i.e., H || c, cannot confirm the 
existence of horizontal node, and clarify the pairing sym- 
metry in heavily hole doped KFe2As2. 

High quality single crystalline samples of KFe2As2 
were grown by a self-flux method which will be de- 
scribed in detail elsewhere. 25 Figure 1 shows the tem- 
perature dependence of susceptibility below 5 K observed 
at fi H = 1 mT after zero filed cooled and field cooled. 
The T c was estimated to be 3.5 K. The ijSK measure- 
ments were conducted both at M15 beamline at TRI- 
UMF, Vancouver, Canada and at RIKEN-RAL Muon 
Facility in the Ruthreford Appleton Laboratory, Didcot, 
UK. The transverse field (TF) ^SR measurements were 
performed at temperatures from 20 mK to 5 K. In two 
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Fig. 1. (color online) The dc magnetic susceptibility at hqH = 
1 mT observed after zero field cooled (ZFC) and field cooled 
(FC). 



sets of TF-/iSR measurements, the magnetic field, was 
applied in parallel and perpendicularly to the crystallo- 
graphic c-axis, respectively, and always perpendicularly 
to the implanted muon spin polarization. The sample 
was field cooled at the measured magnetic fields to min- 
imize disorder of the FLL due to flux pinning. Since the 
muons stop randomly on the length scale of the FLL, the 
muon spin precession signal provides a random sampling 
of the internal field distribution in the FLL state. 

The time evolution of the muon spin polarization sig- 
nals measured with /.loH = 24 mT parallel to c-axis 
above and below T c are presented in Fig. 2(a). Under a 
TF and below T c , implanted muons experience an inho- 
mogeneity of the field due to FLL formation that leads to 
relaxation. As shown in Fig. 2(a), complete depolariza- 
tion is observed at 2 K, indicating that the entire volume 
falls into the superconducting state. Figure 2(b) shows 
the magnetic field distributions P(B) obtained by means 
of the fast Fourier transforms of TF-/1SR signals. For the 
superconductor in the FLL state, P(B) is uniquely deter- 
mined by the magnetic penetration depth A and coher- 
ence length £. 26 The P(B) at 5 K looks very symmetric 
around the frequency /o, which corresponds to the ap- 
plied field, because the applied field penetrates into the 
sample homogeneously. On the other hand, P{B) at 2 K 
becomes asymmetric due to FLL formation. Fig. 2(c) 
shows the calculation results of magnetic field distribu- 
tion for both triangular and square symmetries of FLL. 
The difference between the low field cutoff (I) and saddle- 
point (s) fields is larger in a square FLL than in a trian- 
gular one. Thus, a large shoulder on the low field side of 
the peak is a signature of a square FLL. Comparing the 
observed and calculation results, it is obvious that the 
FLL forms triangular symmetry, consistent with SANS 
results. 22 

The spacial magnetic field distribution in a FLL state 
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Fig. 2. (color online) (a) Muon TF precession signals at normal 
state (5 K) and flux-line lattice (FLL) state (2 K) of single crys- 
talline sample KFe2As2 with fioH = 24 mT parallel to c axis. 
The solid lines represent the result of fitting analysis (see text), 
(b) Local magnetic field distribution P(B) in the normal and 
FLL state. The dotted line stands for the frequency /o, which 
corresponds to the applied field, (c) Theoretical P(B) distribu- 
tions (A = 60 nm, !; v = 10 nm, and fioH = 0.5 T) for triangular 
FLL (green) and square FLL (pink) state, (d) The contour map 
of the FLL at 2 K in real space. 



the cutoff parameter. The TF-^tSR time spectra were 
fitted to a theoretical function by assuming the internal 
field distribution P(B) obtained from Eq. (2) and ac- 
counting for the FLL disorder and the nuclear moment 
contributions, 

A{t) = Ae- {a > a ^l 2 f P(B)e l ^" Bt -^dB, (3) 



where A is the initial asymmetry, a p and cr n denote the 
relaxation due to the FLL disorder and the nuclear mo- 
ment, respectively, 7 M /27r = 135.5 MHz/T is the muon 
gyromagnetic ratio and </> is the initial phase. Fig. 2(d) 
shows the contour map of a FLL at 2 K obtained by the 
fitting analysis. 

Figure 3 shows the temperature dependence of the in- 
verse squared magnetic penetration depth l/A 2 h (T). A 
dip was observed around 0.75 K, similar to the results 
in Bai_ a: K a ;Fe2As2, 13,27,28 which are well described by 
the multiple full-gap behavior. As shown in Fig. 3, both 
single s- and d-wave curves show significant departure 
from the experimental data. The experimental data were 
analyzed within the framework of the phenomenological 
a model by assuming two independent contributions to 
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where K is the reciprocal lattice vector, Bo is the aver- 
age internal field, r is the vector coordinate, and is 



A' 2 (T,A) _ A~ 2 (r,Ai) 
A" 2 (0) ~ W A- 2 (0) 



(1 



w) 



,A - 2 (T,A 2 ) 
A- 2 (0) ! 



(4) 



3 




Temperature [K] HI H 2 



Fig. 3. (color online) Temperature dependence of the inverse 
squared magnetic penetration depth 1/A 2 , obtained with (iqH = 
24 mT parallel to c-axis. The solid black line represents the re- 
sult of a fitting analysis with the two-full-gap model. For com- 
parison, the calculated lines for clean s-wave and d-wecve pairings 
are also presented. 

A- 2 (T,A,) _ 2 r r°° df(E) Esin9dEdO 
A- 2 (0) ~ 1 + J J At (T,9)dE ^£2_ Ai(T; )2 ' 

(5) 

where w is the ratio of gap energy between the two gaps, 
f(E) is the Fermi function, A^(T, 9) is the superconduct- 
ing gap; for simplicity, the Aj(T, 9) = Aj(T) cos 9 for the 
nodal gap and Aj(T) for the full gap. The tempera- 
ture dependence of the gap is approximated by A j (T) = 
A(0)tanh{1.82[1.018(T c /T- l)] 51 }. 30 The result of 
analysis with full gap model is presented in Fig. 3 by 
solid black line, yielding 1/A 2 b (0) = 26.5(1) /Jin~ 2 , w = 
0.75(1), Ai(0) = 0.48(2) meV, A 2 (0) = 0.11(1) meV. 

The field dependence of \ a b observed at T — 0.03T C 
is shown in Fig. 4. It clearly exhibits two kinds of 
field dependence, where the gradient changes at around 
H/H C 2 =0.1. The increase of A a {, is attributed to the 
anisotropic order parameters because it is supposed to be 
a constant in the conventional s-wave superconductors. 
A fitting by the relation X(H/H c2 ) = X(0)[l+rf(H/H c2 )] 
provides a dimensionless parameter 77 which represents 
the strength of the pair breaking effect. From the anal- 
ysis, we find that i] = 1.1(1) and rj = 0.1488(2) for 
H/H c2 < 0.1 and H/H c2 > 0.1, respectively. The large 
value of 77 for H/H C 2 < 0.1 is attributed to be quasi- 
particle excitations from the smaller gap A2 due to the 
effect of Doppler shift. On the other hand, the 77 is small 
above H/H C 2 = 0.1 because the larger gap Ai is robust 
against the field. 

the data with H _!_ c were analyzed using the following 
function, 

A(t) = A s e-^+^ t2 / 2 cos ( ltl Bt + </>) (6) 

where A s is the initial asymmetry and the depolarization 
rates a v is the damping due to the FLL formation. The 
superconducting part of the Gaussian depolarization rate 
a v can be converted into A by 26 

alhl = 0.00371$ 2 A 4 , (7) 

where $0 is the magnetic flux quantum. Figure 5 shows 
A~ 2 (T) obtained from the measured a v (T) by means of 
Eq. (7). It is obvious that T- linear behavior was clearly 



Fig. 4. (color online) Magnetic field dependence of A a j, obtained at 
T = 0.03T C with various applied fields. The solid lines show the 
result of fitting analysis (see text). The H c 2 at this temperature 
is approximately f .8 T. 
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Fig. 5. (color online) Temperature dependence of the inverse 
squared magnetic penetration depth 1/A 2 obtained with fioH = 
16 mT perpendicular to c-axis. The solid and dashed lines rep- 
resent the result of fitting analysis with the line node gap model 
by two-gap and single gap, respectively. 



observed below ~ T c /2, suggesting the existence of line 
node in the superconducting gap. The curve fit assuming 
a T c = 3. IK, which is obtained by the fitting of 1/A 2 h , 
using eq. (4) and (5) with A(T, 9) = A(T) cos6> perfectly 
reproduced the data in Fig. 5 with 1/A 2 = 3.84(8) /im~ 2 , 
w = 0.6(1), Ai(0) = 0.8(2) meV, and A 2 (0) = 0.5(4) 
meV. 

We now discuss possible gap structure in KFe2As2- 
Regarding the results obtained with H || c, both tem- 
perature and field dependence of A indicate the two- 
full-gap or anisotropic gaps. Because the increase of A 
against field is attributed to the anisotropic order pa- 
rameters and the associated nonlinear effect due to the 
Doppler shift of the quasiparticles in the nodal region 
(A(fc) ~ 0). 31 It is predicted that r\ <C 1 for the isotropic 
s-wave pairing because the finite gap prevents the shifted 
levels of quasiparticle excitations from being occupied at 
low temperatures. Therefore, the finite value of r\ means 
the existence of quasiparticle excitations. The field de- 
pendence of A is expected to be stronger when the phase 
space satisfying A(fc) ~ has larger volume, 32 so that 
the quasiparticle excitation observed below H/ H c2 = 0.1 
suggests the collapse of small gap by the field. The 77 be- 
comes small above H/H c2 =0.1, because the quasiparti- 
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cle excitations are suppressed, suggesting the larger gap 
is robust against the field. The value 77 = 1.1(1) is similar 
to the value of r\ = 0.95(1) observed inYNi 2 B 2 C, 33 which 
is a superconducotr having point node in the gap. 34 On 
the other hand, temperature dependence of A~ 2 with 
H _L c shows clear T-linear behavior, suggesting the ex- 
istence of line node in the gap. This is not consistent 
with the results obtained with H \\ c. 

According to the theoretical calculation of the temper- 
ature dependence of the superfluid density n s (1/A 2 oc 
n s ) 7 it is expected that the temperature dependence is 
different from the relation between the direction of line 
node and that of applied field. For example, it behaves 
T-linear and T 3 for the line node perpendicular to the 
field and that parallel to the field, respectively. 35 These 
temperature dependence was observed experimentally in 
B phase of UPt3 , 36 in which temperature dependence of 
A~ 2 parallel and perpendicular to c-axis behaves differ- 
ent. Actually, UPt 3 is expected to have not only line 
node but both line node and point node, but anisotropic 
behavior was observed. Thus, considering the result of 
UPt3, it is obvious that the temperature dependence of 
A~ 2 in KFe2As2 indicate the existence of line node. To 
explain the observed results without any discrepancy, the 
existence of horizontal line node in basal plane is ex- 
pected. This interpretation is perfectly consistent with 
other experimental results, i.e., NQR and specific heat, 18 
thermal conductivity 19 and penetration depth. 21 Re- 
garding the FLL symmetry observed with H || c, the 
interaction between vortices should be isotropic if a line 
node exists in the basal plane. Moreover, the energy scale 
of the obtained small gap ~ 0.11 meV ~ 1.9 T is one or- 
der larger than the applied field of Hj H c2 = 0.1 (H ~ 0.2 
T). The large 77 value obtained below H/H C 2 = 0.1 de- 
spite a small applied field also suppose the existence of 
line node. 

In summary, we have performed the first measure- 
ments of the anisotropic magnetic penetration depth A 
in the iron-based superconductor KFe2As2- The A~ 2 de- 
creases proportionally with T when measured perpendic- 
ular to c-axis but with a higher power of T (~ T 4 ) par- 
allel to c-axis. Considering the theoretical calculation, 
the data can be understood by a superconducting gap 
function with a horizontal line node in the basal plane. 

We thank the staff of TRIUMF's Center for Molecu- 
lar and Materials Science for technical assistance. K.O. 
thanks K. Kubo and H. Sakai for the helpful discussions. 
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